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Possible influence of primordial black hole (PBH) evaporations on cosmic microwave background 
(CMB) is investigated. The spectrum distortions of CMB from the black-body spectrum are de- 
scribed by the chemical potential /i and the Compton parameter y. From COBE/FIRAS limits on 
■ /i and y, the power law index n of primordial density fluctuations and the mass fraction of PBHs (5 

' are constrained by employing the peak theory for the formation process of PBHs. Constraints set 

here are n < 1.304 and n < 1.333 in the thresholds of peaks ("th ~ 0.7 and ("th ~ 1-2, respectively, 
for the PBH mass range between 2.7 x 10"g and 1.6 x 10^^ g, and n < 1.312 and n < 1.343 in 
the thresholds of peaks C^th = 0.7 and ("th = 1-2, respectively, for the PBH mass range between 
1.6 X 10^^ g and 3.5 x 10^^ g, which correspond to the comoving scales between 3 x 10~^* Mpc and 
4 X lO^^'^ Mpc. The constraint on the PBH fraction, which is the direct probe of the amplitude 
of density fluctuations in these scales, stays at almost the same value as /3 < 10~^^ in these mass 
ranges. It is also found that, with these constraints, UV photons injected by PBH evaporations are 
^ , unlikely to ionize the majority of hydrogen atoms. 
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22 ■ I. INTRODUCTION 

^ ' The properties of the primordial density fluctuations in large scales are revealed by Wilkinson Microwave Anisotropy 
^ , Probe (WMAP) and galaxy redshift surveys [H, [1, [I]- The results of these observations are essentially consistent with 
' the predictions of the inflationary model, i.e., primordial density fluctuations with an almost scale- invariant spectrum 
, and the random Gaussian statistics. 
T-H ■ Taking the second order terms of the inflaton potential and density perturbations into account, however, we expect 
] to have the departure from the scale invariant spectrum and the existence of the non-Gaussian components. Therefore, 
. precise measurements of the spectrum and the non-Gaussian components are crucial for a detailed understanding of 
inflation. For example, the first year WMAP data together with galaxy redshift surveys and Ly-a measurements 
OO , preferred the non-zero running spectral index while the third year data with new Ly-a measurements do not 

' strongly favor the running Q. 
^ , The Ly-a measurements are sensitive to the power spectrum at ^ 1 Mpc, while measurements of strong gravitational 
lensing have a sensitivity on much smaller scales such as 10-100 kpc. These methods are, however, limited in both 
scales and epochs. So far the best probe of the small scale density fluctuations is provided by the abundance of 
primordial black holes (PBHs) [Hi ■ The formation of PBHs took place during the radiation dominated epoch due to 
the gravitational collapse of the high density region at the horizon scale when the amplitude of over-density exceeded 
a critical threshold. Therefore the resultant mass spectrum and the abundance of PBHs depend on the amplitude 
of the power spectrum for the primordial density fluctuations at the horizon crossing epoch. It is known that PBHs 
eventually evaporate while emitting Hawking radiation . The lifetime of a PBH is proportional to the cubic of its 
mass. Therefore, PBHs with mass less than lO^^g should have evaporated away by the present epoch. When they 
evaporate, they emit black-body radiation and numerous kinds of particles such as neutrinos, electrons and protons. 

Since PBHs with mass larger than lO^^g survive in the present universe, we can set a constraint on their abundance 
from the fact that PBH density cannot exceed the average matter density observed at present. We can therefore 
introduce (3{M) which is the fraction of the regions of mass M collapsing into PBHs 7]. This can be related to the 
mass fraction of PBHs at the time of formation as P — ppBii.i/Pr,i where ppbha and pr,i are the energy densities of 
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PBHs and background radiation at the PBH formation epoch, respectively. The constraint on the density parameter 
of PBHs at present f7pBH,o < 1 impHes /3 < 10-'^^{M/10^^gy/^ . 

Another important constraint on the abundance of PBHs can be set by the phenomena of evaporation. PBHs 
emit many kinds of particles when they evaporate. Among these are diffuse gamma rays. The emitted photons from 
evaporated PBHs, after the recombination epoch, contribute to the diffuse gamma rays. The upper limit on diffuse 
gamma rays sets a stringent constraint on the abundance of PBHs [1] . It is also known that big bang nucleosynthesis Q 
and entropy production in the early universe also provide constraints on PBH abundance for various mass ranges (see 
the review [1(J|). 

In this paper, we investigate spectrum distortions of cosmic microwave background radiation (CMB) caused by PBH 
evaporation. Naselskii was the first to study the effect of PBHs on the recombination process and to estimate the 
allowed abundance of PBHs that evaporated around the recombination epo ch Jll| . Naselskii and Shevelev estimated 
distortions of CMB due to electrons and positrons from PBH evaporations [l^ • Ricotti et al. have also investigated 
the effects of non-evaporating PBHs on CMB [l^ . X-rays emitted by gas accretion into non-evaporating PBHs were 
found to modify recombination and reionization. Therefore, by studying the effect of such X-rays on the y-distortion 
and the reionization, they obtained a PBH's constraint which is larger than 10^^ Mq. Here, we adopt a modern analysis 
of CMB spectrum distortions, i.e., fi and y distortions on which COBE/FIRAS has set stringent observational limits, 
where /i and y are the chemical potential and the Compton y-parameter, respectively. The process we consider here 
is that photons which have evaporated from PBHs directly induce the /^-distortion in CMB and also hit electrons in 
the surrounding medium. Accordingly, these electrons scatter CMB photons and generate the y-distortion. While 
electrons directly emitted by PBHs also generate the y-distortion via inverse Compton scattering, resultant photons 
are too energetic to be observed as CMB distortions in COBE/FIRAS data. Therefore we ignore electron emission 
from PBHs in this paper. 

Distortions are created when the energy injection into CMB arises between redshift z ^ 10^ and recombination 
z ~ 1000. Before z ^ 10^, CMB can recover the black-body distribution via double Compton scattering or free- free 
emission. After recombination (z < 1000), there is almost no free electron to be a target of photons. Therefore 
we calculate CMB distortions caused by photons emitted by PBHs in this epoch. Consequently, we can obtain the 
constraint on the evaporated PBH abundance using COBE/FIRAS limits, |^| < 9 x 10^^ and y < 1.5 x 10~^ fl^ . 
Note that the mass range of PBHs constrained here is between 10^^ and 10^'^g, since this mass range corresponds to 
the PBHs whose lifetimes are within 1000 < z < 10^. 

In order to estimate the number density of PBHs, we need the threshold value of density fluctuations required for 
the overdense regions to collapse into PBHs. The threshold value was traditionally calculated as the critical density 
contrast in the simplified cosmological model by the analytic method, where the value obtained is (5c = 1/3 for the 
radiation dominant era [5] . There has been resent progress in estimating the threshold value by employing numerical 
simulations based on general relativity [Tsl. [T6|. In this paper, therefore, we adopt the threshold value of Shibata and 
Sasaki 

This paper is organized as follows. In Sec. II we review the formula for evaluating the PBH mass function by using 
the peak theory. In Sec. HI we calculate CMB distortions, fi and y, caused by Hawking radiation from PBHs and we 
set constraints on the primordial power spectral index and the PBH abundance. Conclusions are presented in Sec. 
IV. In this paper we assume that h = 0.70 {Hq ^ h x 100 km/s ■ Mpc), fl^h'^ = 0.022 and fluh'^ = 0.11 Q. And we 
set c ^ h = 1, where and c and h are the speed of light and Planck's constant over 27r, respectively. 



II. PBH MASS FUNCTION 



In this section, we calculate the distribution of PBHs produced from the primordial density fluctuations by following 
the procedure of Green et al.[l3. We adopt the threshold value of density fluctuations required for the overdense 
regions to collapse into PBHs in the horizon crossing epoch obtained by Shibata and Sasaki [l^ in which they 
performed numerical simulation of evolution for the metric perturbations on the uniform-expansion hyper surface. 
Green et al. [TtI related this threshold to the gauge-invariant curvature perturbation on the uniform-density hyper 
surface, (, which is defined as [l3| 

C-7^-i/^, (1) 
p 

where p is the background density, Sp is the perturbed density, TZ is the curvature perturbation and dot represents the 
time derivative. Shibata and Sasaki's threshold depends on the environment of the PBH formation. The threshold 
value ^th is Cth — 0.7 when the density peak is surrounded by a low density region. In contrast, if the peak is 
surrounded by a flat Friedman-Robertson- Walker region, Cth = 1.2. 
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Green et al. advocated the peak theory for the calculation of the PBH mass function because Shibata and Sasaki's 
result is described as the constraint on the peak value of fluctuations. We assume a power-law primordial power 
spectrum, Vn = TlAk / koY'^^ . From WMAP and galaxy survey results, Tic is (24.0 ±1.2) x 10" i° at the scale 
fco = 0.002Mpc~^ We smooth the density fields by a Gaussian window function with comoving size R. Under 
these assumptions, the peak theory gives the comoving number density of the peaks which are higher than v as [l7| 



n{v, R) 



1 (n- 1)3/2 



(z.2-i)exp(-^ ) , 



(27r)2 63/2i?3 

where we employ the high peak limit ^ 1. We relate v to the threshold value of the PBH formation by 



2(fcoi?) 



7^e^((n - l)/2) 



1/2 



th- 



(2) 



(3) 



We can assume Eq. ([2|) to be the comoving number density of PBHs which arc formed from the overdense regions 
with scale R. 

Let us relate the smoothing scale R to the PBH mass. The PBH mass depends on the initial environment around 
the peak as well as the threshold value. For simplicity, however, we assume that PBHs with the horizon mass are 
produced when the overdense regions enter the horizon. The epoch when scale R crosses the horizon is evaluated as 



(4) 



From the comoving entropy conservation, we obtain the relation between the temperature T and the scale factor a 
(which is normalized at the present epoch as oq = 1) 



gV^aT = const. 



(5) 



where 17* is the number of relativistic degrees of freedom. In the radiation era, the horizon mass at horizon crossing 
of the comoving scale R is described from Eqs. ^ and ([5]) as 



MbhIR) 



4n /SttG 
"3 



Hq^Im f g*cq 

1 + ^cq V g* 



1/3 



1/2 



R' 



-1/6 



R 



6.2 X 10s cm 



(6) 



Here, we assume that the redshift of the matter-radiation equality is Zoq — 3200 [ij and the number of relativistic 
degrees of freedom at the equality is g^oq = 3.36. 

The formation epoch of a PBH with the mass A/bh, which we assume as the horizon crossing epoch of this mass 
scale, can be described in terms of the scale factor as 



aBH(A^BH) 



^) Hon]/' 



BH 



X 



1/2 



(7) 
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47r /SttG 
"3 



1 + Zca 



g^c(\ 
g* 



1/3 



1/2 



(8) 



Moreover, employing Eq. ([5]), we can rewrite Eq. (21) as a function of the PBH mass n(z/, M). 

To describe the PBH abundance, one often uses the fraction of regions of mass M which collapse to PBHs, (3{M). 
Obviously, (3{M) also represents the mass fraction of PBHs at the time of formation. The energy density of PBHs 
with mass Mbh at horizon crossing of scale R is expressed as pbh(-^^bh) = -^^BHf^('^, -^^BH)aBH(AfBH)'^- The mass 
fraction of PBHs (3 can be written as 



/3(M] 



BH) 



/Obh(Mbh) 

p 



1 



-1) 

63/2 



3/2 



-{v"^ - 1) exp 



(9) 



Here we adopt the high peak approximation, Eq. 
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One could conclude that the total number density of PBHs with mass Mbh at a given epoch can be described as 
n(i/, Mbh)- However, this is the number density of PBHs when they were formed. For the total number at a given 
epoch, which we express as riBnii^, Mbk, z), we need to subtract the number of PBHs which have been absorbed 
by larger PBHs in subsequent epochs. To take into account this effect, we consider the continuity equation for the 
PBH number density. The comoving number density of PBHs at z, in the range from Mbh to Mbh + c^Mbh can be 
described, in the high peak limit, as 



nBH(l', MBH)rfMBH 



dn{iy, Mbh) 



dM, 



BH 



rfMBH 

1 /'X{n-1)Y^'^ 



47r2A/BH V 6AfBH 



^^^-^i/4(n, Mbh) -^in- 3)i/2(n, Mbh) - 3 



exp 



z/2(n,A/BH) 



^^y-^M^) ^-(-^^^BH)exp(^ j 



(10) 



BH) 



2 (A;2Mbh/^ 



(n-l)/2 



7^,^((n-l)/2) 



1/2 







th- 



fll) 



At the third step in Eq. PH)) . we take the high peak limit, ^ 1. 

Let us discuss the epoch of first PBH formation. For simplicity, we set the end of inflation at this epoch. Accordingly, 
the initial condition of PBH formation can be characterized by the reheating temperature. Under this assumption, 
the minimum mass of PBHs depends on this temperature. The comoving horizon scale at this epoch is expressed 
from Eqs. ([4]) and ^ as 



1 

oH 



g* 

g*eq 



1/3 



1/2 



J- cq 



= 2.0 X 10' 



viooy Vio**Gevy 

where T^h is the reheating temperature. By using Eq. ([6]), we obtain the minimum mass of PBHs as 

AfBH,min(i;,0 - 9.8 X 10^^ (y^^^) §• 



(12) 



(13) 



In our assumption, therefore, n{i', M) = if M < MBH,inin(7rh). 



III. CONSTRAINTS FROM THE CMB DISTORTIONS 



PBHs inject energy into CMB through Hawking radiation. In the early universe, even when there is a large amount 
of energy injection, CMB achieves black-body spectrum by photon-electron interaction, i.e., Compton and double 
Compton scatterings [l^. However, decouplings of these interactions occur one by one after z ^ 10^ so that the 
distortions from the black-body are produced via the energy injection. 

First, the double Compton scattering is decoupled at z ^ 10^. After this decoupling, the total photon number 
between CMB and the injected photons conserves. Accordingly, CMB can no longer achieve the black-body spectrum 
if there are any photon injections after the decoupling of the double Compton scattering, although the injected photons 
are still thermalized by the Compton scattering. The CMB spectrum in this thermal equilibrium state is described 
as the Bose-Einstein spectrum with the chemical potential fi which stands for "distortions from the black-body 
spectrum" . 

Following the double Compton decouphng, the decoupling of the thermahzation comes at z ^ 10^. Due to the 
expansion of the universe, Compton scattering does not work effectively and can no longer establish the thermahzation 
when the time scale of the thermahzation becomes longer than the Hubble time. We must describe the distortion 
from the black-body spectrum due to the photon injections as the Compton-y parameter once the decoupling of the 
thermahzation takes place. 
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These distortion parameters are constrained by COBE/FIRAS, < 9 x 10~^ and y < 1.5 x 10~^. In this section 
we discussion constraints on the PBH abundance or the spectral index of density fluctuations by the limits of the 
chemical potential and the Compton-y parameter. 



A. ^-distortion 



When energy Q is injected continuously into CMB, time evolution of the chemical potential, which describes the 
distortion of the CMB energy spectrum, is given as [20| 



dt t-Dc{z) 



Q 

1.4—. 

Pi 



Here ^dc is the time scale for the double Compton scattering 

ioc 2.06 X 10^3(1 - Yp/2)-^{n^h^)-^z-'^/'^ s, 
where Yp is the primordial helium mass fraction. The solution of Eq. ([T4|) is given by Hu and Silk 20] as 



/i = 1.4 / at — exp 

Jo Pi 



z 

ZDC 



(14) 



(15) 



(16) 



ZDc = 1-97 X 10" 1 - - 



1 / 



2 V0.24 



-2/5 



0.0224 



-2/5 



(17) 



^freeze = 2.86 X 10" ( 1 - - 



1 / 



0.24 



-1/2 



0.0224 / 



-1/2 



(18) 

reczo is the redshift 



where zdc is the characteristic redshift for decoupling of the double Compton scattering and Zf, 
when the injected energy can no longer be thermalized. 

As mentioned before, the /i-distortion takes place if there are any photon injections during the epoch after decoupling 
of the double Compton scattering and before decoupling of the Compton scattering, i.e., Zfreezc < z < zdc- 

Let us now estimate the energy injection Q due to the Hawking evaporation of PBHs. A Schwarzschild black hole 
with mass M emits particles with spin s and total energy between E and E -I- dE at a rate per degree of freedom 



dE = 

dtdE 2TTh 



-1 



i2s 



dE, 



(19) 



where T{M) is the temperature of a black hole with mass M and written as 

1 



r(M) 



SttGM 
f M 



1.0 



Vio 



13n 



GeV. 



(20) 



Here T ^ is the dimensionless absorption probability of emitted species. For simplicity, we assume that only photons 
contribute to Q. In the case of photons, is described as 



64G"*M4£;4/3, E < fcT(m), 
27G^]VPE^, E > fcT(m). 



(21) 



The PBH mass AI decreases due to the Hawking evaporation. The mass loss rate of a PBH can be written as [22 



dM Tt; . ^ T 1 

— = -5.34 X 10^^f{M)M-^g sec"\ 



(22) 



6 



where f{M) is a function of the number of species which are directly emitted, and can be described by the fitting 
formula as 



/(M) = 1.569 + 0.569 



, 0.0234\ ^ /-0.066\ ^ /-0.11\ /-0.394 
'^'^P rn/ + D cxp + 3 exp + exp I 



-0.413 



mv ; + ^ 1, Y{M) ) ^ '''''' [jW) J [nm 

o f -22 \1 ^^^^ /-0.10\ 
TTTT^ + 3 exn ^.TT^ + 0.963 exn .rrr^ 



Here T(M) is written in the unit of GeV as in Eq. (l20l) . 

Since f{M) is a weak function of M, we can approximately integrate Eq. and obtain the time evolution of the 
PBH mass with an initial mass Mbh by 

M(Mbh, t) « [Af-^H - 1-5 X 102V(M)t] , (24) 

where the typical values of f{M) are 1.0, 1.6, 9.8 or 13.6 for M > IQi'^g, M = lO^^g, lO^^g or 10"g, respectively. 
From Eq. ([22]). we can also obtain the life time t of the PBH with the initial mass Mbh , 

/A'/bh A/*^ 
'^^^7(M) ^^^^ 



Finally we can write the energy injection rate due to PBHs as 

fM„{t) poa 11^ 

'M„i„(t) 



j-M„(t) poo ij^ 

Q{t)^ dM^a dEa-'{t)njiniMBH)-j^{M{MBH,t))E, (26) 



where Muit) is the horizon mass at t and M-aiin{t) is the minimum initial mass of the PBHs with a lifetime t = t. 
The M„iin(t) can be evaluated by taking the inversion of Eq. 

Now we are ready to calculate the evolution of the chemical potential induced by photon injections from PBHs. 
Let us first assume that the lifetime of the smallest PBHs which are formed right after the inflation epoch is shorter 
than the time scale of the double Compton decoupling ^dc • Under this assumption, we expect the existence of PBHs 
that evaporated away at the epoch of the double Compton decoupling. The mass of these PBHs is obtained as 
Mdc = 2.7 X lO^^g, which corresponds to the mass scale of the horizon when the temperature Tdc = 3.1 x lO^^GeV. 
If the reheating temperature of inflation is higher than this temperature, i.e., Tih S> Tbc the mass of the smallest 
PBHs formed after inflation is less than Mdc- 

Under this assumption, we can compute the chemical potential ji by using Eqs. (fTi|) and ([^5]) . We show the 
resultant chemical potential as a function of the primordial spectral index of the power spectrum in FIG. [1] The left 
and right panels correspond to the cases of the critical threshold of PBH formation Cth = 0.7 and C,t\\ — 1-2. The 
COBE/FIRAS upper limit on the chemical potential is /i < 9.0 x 10^^. We acquire the constraints on the spectral 
index n as n < 1.304 for C,th — 0.7 and as n < 1.333 for Cth = 1-2. 

Fig. [2] shows the production rate of the chemical potential per PBH mass. The thick solid line is the chemical 
potential per PBH mass for Cth = 0.7 and n — 1.304 and the thick dotted line is for Cth — 1-2 and n — 1.333. The 
vertical thin dashed line is for A/dc- We can see that the majority of contributions to the chemical potential is 
made by PBHs with mass around Mdc- This may look counterintuitive because Eq. p4p implies that, the larger 
the PBH mass is, the larger the chemical potential per evaporated PBH is produced. However the number density 
of PBHs with a smaller mass is much larger than that of PBHs with a larger mass. From the comparison of these 
two contributions, we find that PBHs with the smallest mass existing at zdc, which is Mdc: provide most of the 
contributions for the power spectrum with n > 1. Accordingly, PBHs with masses below Mdcj which evaporated 
away before zdC; are not relevant for the /i-distortion as is shown in Fig. [2l Since the reheating temperature only 
determines the minimum mass of PBHs, the constraint on the spectral index from the /Lt-distortion does not depend 
on the value of the reheating temperature as far as Tjh 3> Tbc- 

The other characteristic PBH mass scale is Mftcczc = 1-6 x lO^^g, which corresponds to the mass scale of the 
horizon when the temperature Tfrcozo = 7.6 x lO^^GeV. PBHs with Affreczc evaporate away at the epoch of the 
thermalization decoupling Zfroczc- This mass scale is represented as a vertical thin solid line in Fig. O We can see 
that the contribution of PBHs on the chemical potential /i is quickly suppressed once the mass exceeds Mfroozo- This 
feature, which is represented by the kink in each thick line in Fig. [5] at Affioezc, can be easily explained by the fact 
that PBHs with mass below Mfreeze can only partially evaporate by Zfroczc, and contribute little to the /i-distortion. 

Let us next consider the case of the reheating temperature Tj-h ^ Tbc- In this case, the mass of the smallest 
PBHs formed after inflation, which is described as MBH,min, is greater than Mdc- Accordingly, there is almost no 
evaporation of PBHs at zdc- Instead, the evaporation takes place at a later epoch. Generally speaking, the lower 
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1.30 1.302 1.304 1.306 1.308 1.33 1.332 1.334 1.336 



spectrum index n spectrum index n 

FIG. 1: The chemical potential for each spectral index. We assume the reheating temperature Trh ^ Tuc (see text). The 
critical thresholds of PBH formation are taken as C^th = 0.7 and 1.2 in the left and right panels, respectively. The left and 
right panels correspond to the cases of PBH formation at the density peaks surrounded by the average and low density regions, 
respectively. The dotted line is the upper limit of the COBE/FIRAS observation, /i — 9.0 x 10~^. The region under the dotted 
line is allowed. We find the upper limit on the spectral index n as n < 1.304 for ("th ~ 0.7 and as n < 1.333 for ("th ~ 1.2. 



the reheating temperature is, the later the PBH evaporation takes place. Therefore, we expect that the constraint 
on the power law index n becomes looser for a lower reheating temperature. This behavior is shown in Fig. [S] Note 
that if the reheating temperature is so small that MBH,min(T'rh) > .^freeze, the majority of PBHs evaporate away after 
^freeze. Accordingly, the constraint on the power law index becomes much looser as is shown in the left hand side of 
the vertical soHd Kne of Fig. O 

Now, we will constrain the mass fraction of PBHs (3. If T^h ^ Tbcj wc find that greatest contribution to /i is made 
by PBHs with the mass AfBH.min. In this case, therefore, the constraint on the spectral index for a given reheating 
temperature shown in Fig. [3] can be converted into the mass fraction of the PBHs with the minimum mass. Describing 
the reheating temperature by AfBH,min with Eq. ([TH)) . we obtain the constraint as shown in Fig.[31 On the other hand, 
if Trh 3> Tuc, the most efficient PBHs for the ^-distortion are those with a mass of A/dc- Therefore the constraint 
from /i-distortion does not provide sufficient clues for understanding the mass fraction on the scale below Afnc- In 
Fig. [21 however, we extrapolate our procedure to slightly smaller masses, since contributions from masses slightly 
smaller than Afoc may still be important as is shown in Fig. [51 The constraint we obtain is: /? < 10~^^ between 
A^DC < Af < Afficeze, which is tighter than previous constraints found in (lo| . 




log[MBH[g]] 



FIG. 2: The production rate of the chemical potential per PBH mass. We assume the reheating temperature Trh 'S> Tuc (see 
text). The thick solid and thick dotted lines correspond to the cases with i^th ~ 0.7 and n — 1.304, and ^th ~ 1.2 and n = 1.333, 
respectively. The vertical thin dashed line represents the characteristic PBH mass scale Muc while the thin solid line represents 
A/frcczc. It is clear that PBHs with the mass around Afoc give a dominant contribution. 



B. Compton y-parameter 

After Zfreeze, the injected photons are no longer thermalized by the Compton scattering. These photons scat- 
ter electrons and the resultant energetic electrons scatter the CMB photons through the inverse Compton process. 
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FIG. 3: The constraint on the spectral index as a function of the reheating temperature from the /i-distortion. The thick solid 
line and the thick dotted line are the upper limits with the PBH critical threshold ("th ~ 0.7 and ("th = 1.2, respectively. The 
vertical thin dashed and thin solid lines represent the reheating temperatures whose horizon scales correspond to PBHs with 
masses Muc and A/frcozo, respectively. 




FIG. 4: The constraint on the PBH abundance from the /i-distortion. The allowed region is under the thick line. The vertical 
thin dashed and thin solid lines represent the locations of Muc and Af freeze, respectively. 



Accordingly, Compton y-distortions are induced. 

When Te T, where Tg and T are electron and CMB temperatures, respectively, the Compton y-parameter is 
described from the Kompaneets equation as [1^ 

y = [ dt^-^Uear, (27) 

where rig and ax are the electron number density and the Thomson scattering cross section, and ^(zfrcozo) Etnd tQ are 
the freeze out time and the present time, respectively. 

The electron temperature Tg is controlled by Compton heating and Compton cooling as 



t—kTeUe = near E^n-^ - Ane(JTk{Te - T)— [kTj I — I - 15-t 

OiT ^TT'C -LO \ TTIq I Qj 



-df = me ""^^^"-^^ + E,n,n.ar + -E,— - 2-n, + Sn,{E,), (29) 



where n^{E^) is the injected photon number density with the energy being from E^^ to E^ + dE-^, and 5n^f{E^) is the 
production rate of injected photons. The first term of the right hand side (rhs) of Eq. describes the Compton 
heating of electrons by injected photons, the second term describes the Compton cooling of electrons from CMB 
photons and the third term describes the adiabatic cooling due to cosmic expansion. The first two terms of the rhs 
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of Eq. (|29|) describe the Compton cooling for injected photons, and the third and forth terms describe the adiabatic 
coohng. 

For obtaining the electron temperature Te, we employ the analytic approach in Ref. (23 |. When we assume a high 
electron temperature (Tg ^ T) and the steady state {dT^/dt = 0), we obtain the electron temperature from Eq. 
as 

m = ^^^^^E,n,. (30) 

Here we ignore cosmic expansion. 

Integrating Eq. (|29p over the photon energy, we can obtain the equation for the total photon energy density, 

d /"^ ,„ „ /"^ ,„ E^-AkT, 

dt 



dEjEjUj — ~ / dEj — EjU^neCTT + Q, (31) 



where we ignore cosmic expansion again because the time scale of the Thomson scattering is shorter than the cosmo- 
logical time. The second term Q = J dE^E^ {5n^{E^) / dE^) on rhs is the total injected energy rate and is described 
as Eq. (|26p because we consider the PBHs as the only source of the energy injection. The first term on the rhs of Eq. 
pip describes the energy loss rate per unit time and, hereafter, we express this term as E\oss- From Eq. (j30p . we can 
rewrite the electron temperature in terms of £^ioss as 

fcTe = — (32) 

ne<JT 

In order to acquire the electron temperature, we need to estimate the energy loss rate -Eioss- A photon with energy 
E^ loses the energy E^{E^ — 4:kTe)/me per Compton scattering. Hence, the energy loss of a photon with high initial 
energy E-yQ ^ kT^ within Hubble time is approximated as 

SE,^^"^. (33) 

Substituting the Hubble time and the typical photon energy from PBHs, our present target of interested, we obtain 
SE-y ^ £'70- This implies that the injected photon energy is fully transferred into electrons in the Hubble time. 
Therefore we can approximate the energy loss rate -Eioss as i?ioss — Q- Substituting Q for i?ioss in Eq. (|32p . and Eq. 
(|32ll for Eq. we obtain 

y= / dt^, (34) 

where the upper bound of the integration t{zj:cc) is the recombination epoch, which is introduced since the injected 
energy does not transfer into the background electrons once the optical depth becomes very low after recombination. 

We calculate the y-distortion for each primordial power spectral index under the assumption that the reheating 
temperature is much higher than Tfj-eczc = 7.6 x lO^^GeV which corresponds to Affi-ooze- The result is shown in Fig. 
[SJ The upper limit on the y-parameter obtained by COBE/FIRAS is ?/ = 1.5 x 10~^ which is the dotted line in Fig. 
[5j Therefore, we obtain the constraint on the spectral index n < 1.312 for Cth = 0.7 or n < 1.343 for Cth — 1-2. The 
y-parameter per PBH mass is plotted in Fig. [HI The thick solid line is the y-parameter per PBH mass with Cth — 0.7 
and n — 1.312, and the thick dotted line is that with Cth — 1-2 and n ~ 1.343. The vertical thin solid line is the 
location of Mfreozc and the vertical thin dotted line is Afnc = 3.5 x 10^'^g. The PBHs with Mrc evaporate away at 
the recombination epoch, z w 1000. Because PBHs with masses smaller than Affroczc have evaporated away before the 
redshift Zheeze, they provide no contribution to the y-distortion. It is shown in Fig. [5] that the largest contribution to 
the y-distortion is from PBHs with mass Mfrcczc- 

Next, we calculate the y-distortion in the case of the reheating temperature Tih ^ Tfrcczc- We plot constraints on 
the power law index n for a given reheating temperature in Fig. [T] The thick solid and thick dotted lines are the 
constraint of the spectral index for Cth = 0.7 and Cth — 1-2, respectively. The vertical thin solid and thin dotted lines 
represent Tfrcozo and Trc = 1-7 x lO^GeV which corresponds to M^c, respectively. When the reheating temperature 
is smaller than Tfreeze, the larger spectral index is allowed because contributions on the y-distortion from PBHs with 
masses between AfBH.minCTrh) and Af freeze are missing similar to the case of the chemical potential. 

We calculate the upper bound of p from the constraint on the spectral index which is plotted in Fig. [T] We plot 
the result as the solid line in Fig. [51 As a reference, we plot the constraint from the chemical potential as the thick 
dotted line. These constraints, from /^-distortion and y-distortion, are complementary, as is shown in the figure. On 
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1.308 1.31 1.312 1.314 1.316 1 338 1.34 1.342 1.344 1.346 1.348 

spectrum index n spectrum index n 

FIG. 5: The y-distortion as a function of the spectral index. We assume Trh S> Tfrcozo- The critical thresholds of PBH formation 
are taken as i^th ~ 0.7 and 1.2 in the left and right panels, respectively. The dotted line is the upper limit of the COBE/FIRAS 
observation. The region under the dotted line is allowed. The upper limits on the spectral index n are n < 1.312 for ("th ~ 0.7 
and n < 1.343 for Cth = 1-2. 



the mass scales Mfreozo < -^bh, the y-distortion provides a stringent constraint, while the /^-distortion gives tighter 
constraint on Mfreeze > Mbh- We find that the upper bound of /3 is 10~^^ between Mfreeze < M < Mrc- 

Let us explain the above constraint on /3 from the y-distortion more intuitively. Most of the j/-distortion is produced 
by PBHs with minimum mass at the redshift Zftoczc- Here we denote this minimum mass as Mmini. If T^h ^ Tftc 
A^mini = A%oczc whilc Mniini ^ AffrcozG for T^h < Tfreozc- Wc cau approximate Eq. as 



iccze; 



,.9iM^, ,35, 

where Q(Mmini) is the total injected energy between Zfrcozc and Zc, which describes the epoch when the PBHs 
with Mmini are evaporated away. The maximum injected energy from the PBHs with Afmini is Ppbh (.A/mini). The 
temperature at formation of PBHs with A/bh is Tbh ~ 3 x 10^GeV(A/BH/10^^)^''^ from Eq. p2|) so that the formation 
redshift leads to zbh ^ 10^^(A/bh/10^'^)^/^. Since /3 is defined as the fraction of PBHs at the formation epoch, we 
must take into account the evolution due to cosmic expansion to estimate the PBH fraction at Zg. PBHs can be 
treated as a matter component while the universe is dominated by radiation. Therefore the evolution factor between 
zbh and Zc can be written as (zbh + l)/(-^c + 1)- Eventually we can relate the y-parameter to P as 

^ PPBH (A/mini) f^PBH ^bh + 1 ^ /o^\ 

where the inequality is introduced by the fact that the maximum injected energy from PBHs provides only the upper 
bound for the y-distortion because PBHs can produce not only photons, but many kinds of particles, i.e., electrons, 
neutrinos, and so on. This relation is drawn in Fig. [8] as a thick dashed line. We find this intuitive approach is 
consistent with our previous constraint. 



IV. CONSTRAINT FROM THE OPTICAL DEPTH 



In this section, we investigate the possible reionization of neutral hydrogens after recombination due to evaporation 
of PBHs. 

We adopt the cumulative number of ionizing photons which arc produced by PBH evaporations as the criterion 
of the ionization. For simplicity, we ignore the cascade decays of high energy photons which can result in a larger 
number of photons. If we take this effect into account, the constraint on /3 which we will obtain below, shall be much 
tighter. 

To achieve a volume- weighted ionization fraction of 99 percent, it is found that 5 to 20 cumulative photons per 
hydrogen atom are needed 25, 26]. The particle emission rate of PBHs is given by Eq. (fT9|l so that the total 
cumulative ionizing photon number per hydrogen atom for given time t is described as 

1 l-MHit) j-OO 

= — / dMnu dE nBHiMBH)^^{M{MBn,t)), (37) 

na TT-H JA/„i„(t) Jo dtdE 

where riu is the comoving hydrogen number density. Adopting the upper bound of the power law index n from 
the previous section as n ;< 1.3, we find that the cumulative ionization photon number is so negligible ('~ 10~^ per 
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FIG. 6: The production rate of the j/-parameter per PBH mass. We assume Trh 3> Tfroczc. The thick soUd and thick dotted 
hues correspond to the cases with ("th = 0.7 and n = 1.312, and ("th = 1.2 and n = 1.343, respectively. The vertical thin solid 
and dotted lines represent Mfreeze and Mrc, respectively. It is shown that the contribution of the PBHs with mass Mf,cczc 
gives a dominant contribution so that the abundance of the PBHs with the mass Mfrcozc- 
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FIG. 7: The constraint on the spectral index as a function of the reheating temperature from the j/-distortion. The thick solid 
line and the thick dotted line are the upper limits with the PBH critical threshold ^th ~ 0.7 and ^th = 1.2, respectively. The 
vertical thin dotted and thin solid lines represent the reheating temperatures whose horizon scales correspond to PBHs with 
masses Mrc and Mfreeze, respectively. 
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FIG. 8: The constraint on the PBH abundance from the j/-distortion (thick solid line). The allowed region is under the thick 
solid line. The vertical thin solid and thin dotted lines represent the locations of A/froozc and Afac, respectively. The thick 
dotted line is the constraint on the PBH abundance from the /i-distortion. The thick dashed line is obtained by Eq. (|36p . Here 
we adopt the COBE/FIRAS upper limit for y, and describe Zc and zbu in terms of A/bh. 



12 



hydrogen atom) that reionization cannot be caused by ionizing photons injected by PBHs unless the cascade decays 
of high energy photons take place very efhciently such that one energetic photon decay into, say, 10^ ionizing photons. 

V. CONCLUSION 

In this paper, we investigated the possible influence of PBH evaporation on CMB. We set new and stringent 
constraints on the PBH mass fraction and the primordial power spectrum of density fluctuations from the observational 
upper bounds of /i and y-distortions obtained by COBE/FIRAS. For this analysis, we employed the peak theory, with 
a threshold described as Cth- 

It was shown that fj, and y-distortions set limits on different ranges of PBH masses. From the /i-distortion, we can 
test the mass range between 2.7 x lO^^g and 1.6 x lO^^g, which correspond to Mdc and Mfroozo, that is, the masses 
of PBHs which evaporate away at the decoupling epochs of the double Compton and the thermalization, respectively. 
We obtain constraints on the power law spectral index as rt < 1.304 and n < 1.333 for Cth — 0.7 and Cth = 1-2, 
respectively. For the PBH abundance, we set the limit /3 < 10~^^ for the mass range we considered here. 

On the other hand, from the y-distortion, we investigated the PBH mass range between 1.6 x lO^^g and 3.5 x lO^^g, 
which correspond to Afficozo and M^c- Here Mrc corresponds to the mass of PBHs which evaporate away at the 
recombination epoch. We obtained constraints on the power law spectral index as n < 1.312 and n < 1.343 for 
Cth = 0.7 and Cth = 1.2, respectively. For the PBH abundance, we set the limit /3 < 10~^^ for the mass range 
considered here. 

It turns out that our constraints on the spectral index are looser than those obtained from the observation of CMB 
temperature fluctuations by WMAP satellite and other large scale observations such as galaxy surveys of 2dF and 
SDSS, and surveys of Ly-a forest, which imply a spectral index n = 0.947 ± 0.015 [l|. However, it is suggested 
by many authors that general inflation models produce not a simple power law spectrum but a spectrum with a 
running spectral index which has many branches. Accordingly it may be no use to constrain the power law index from 
observations on various scales. All we can do is to set constraints on the amplitude of each scale associated with the 
observation. From this point of view, the mass fraction of PBHs /3 can provide unique information for the fluctuation 
amplitude on very small scales, Mdc < A/ < M^c which corresponds to the comoving scales between 3 x lO^^^Mpc 
and 4 x 10~^^Mpc. Since these small scales correspond to the Horizon scales right after the end of inflation, we can 
say that the last stage of inflation can by revealed by PBHs. 

Finally, we would like to mention that PBHs can be formed not only by the primordial density perturbations but 
also by the collisions of bubbles of the broken symmetry phase [l^ or by the collapse of cosmic strings [28]. Our 
constraint on /3 is even applicable for such PBHs. 
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